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Abstract 

We investigate the equivalence between spectral characteristics of the Laplace op- 
erator on a metric graph, and the associated unitary scattering operator. We prove 
that the statistics of level spacings, and moments of observations in the eigenbases 
coincide in the limit that all bond lengths approach a positive constant value. 



1 Introduction 

Quantum graphs have attracted much attention in recent years due both to their ap- 
plicability as physical models, and their interesting mathematical properties. We refer 
the reader to |1] and the forthcoming volume ^ for a pot-pourri of new results. Recent 
reviews, dedicated to quantum graphs, include [3l[l]. 

In this article we focus on one feature of interest which is the use of graph models to 
probe the the universality of quantum systems. One of the unsolved paradoxes of quantum 
mechanics is the observation that a great many quantum systems are remarkably similar 
when one makes statistical observations in the semi-classical regime. This manifests itself 
both in the energy levels, and associated energy eigenfunctions. Despite a great deal of 
effort, this universality is poorly understood mathematically. Generic quantum graphs 
exhibit this universal behavior, and represent the most likely system for which a full 
mathematically rigorous proof of this universality will first be found. Important steps in 
this direction have been taken in [5]. 

A quantum graph can be defined in two different, but related, ways (a complete de- 
scription appears in the following section). One may consider a self-adjoint realisation 
of the Laplace operator, or a scattering matrix approach. Mathematically, the scattering 
approach appears to be more tractable, and has formed the basis of most investigations 
[SI El [SI [9l [To] . Moreover, the spectra arising from these quantizations are subtly different. 
This is not as confusing as it might seem, since the statistical properties of both versions 
of the spectrum are believed to coincide when averaged over a large interval. It is the 
purpose of this article to put a concrete mathematical foundation behind this belief. 

The plan of the article is as follows: In the next section we give precise definitions of 
the two ways to describe quantum graphs, and in section [3] describe our main results. In 
section [4] we describe the tools used and then present in section [5] the proofs of our results. 



2 Two descriptions of quantum graphs 



For both constructions of a quantum graph we begin with a graph G = (V, B) where V 
is a finite set of vertices (sometimes referred to as nodes), and B is the set of bonds (or 
edges). Each bond b has a positive length, denoted Lb- The total number of bonds is B. 
We denote by the degree of the vertex f € V, which is the number of bonds emanating 
from it. 



2.1 The Laplace operator approach 

The first way to define a quantum graph is to identify each bond b with the interval [0, L^] 
of the real line and thus define the L^-space of functions on the graph. Then one can 
consider the eigenproblem 

d^ 

- '^^b{x) = X^Ubix). (1) 

This setup has a long history of being used in physical models [HI [12l [13] . It was studied 
by mathematicians since at least the 1980s [H [El [IQ [13 [lEl [H [20]. The study of 
spectral statistics of quantum graphs was initiated in |2H I22j. 

To make the operator in ([1]) self-adjoint one needs to impose matching conditions on 
the behavior of u at the vertices of the graph. One possibility is to impose Kirchhoff 
conditions0 we require that u is continuous at the vertices, and that the probability 
current is conserved, i.e. 

V = for all V eV, (2) 

dx 

v&b 

where the sum is over all bonds that originate from the vertex v and the derivatives are 
taken at the vertex v in the outward direction. The admissible matching conditions were 
classified in, among other sources, [231 [231 [25] . 

A solution to the eigenvalue equation ([2|) on the bond b, can be written as a linear 
combination of plane waves, 

Mxb) = Cfee''^^" + CfeC-'^^" . (3) 

A solution on the whole graph can be uniquely defined by specifying the corresponding vec- 
tor of coefficients c = (ci, . . . , cb,ci, . . . , cb)'^- The elements of the vector c are naturally 
associated with directed bonds of the graph G. 

Imposing the matching conditions, we find that is an eigenvalue if and only if 
S{X)c = c for some explicitly given matrix ^(A). For instance, in the case of Kirchhoff 
conditions given above, 

5(A) = e'^^So, (4) 
where the elements of the 2B x 2B matrix 5*0 are given by 

\"'V2 



^sometimes called "Neumann" conditions 



2 



for {vi, V2) and (fs, ^4) directed bonds of the graph. By L we denoted the vector of bond 
lengths L = (Li, . . . ,Lb) and defined, in a shght abuse of notation, a 2B x 2B diagonal 
matrix 



\ 



(5) 



where x S M^. Note the doubling of dimension which signifies the move from undirected 
to directed bonds. 

The general conditions for factorization ^ are given in \26\ I27j . Even when these 
conditions are not satisfied and S depends on A in a non-trivial fashion, the above fac- 
torization is valid in the limit X ^ 00. Thus, for the purpose of studying the spectral 
statistics it is not a strong restriction to assume that ([4]) applies, as we will do below. 



2.2 The scattering approach 

The second construction considers wave propagation on the graph where each vertex is 
treated as a scatterer and propagation along the bonds is free. This construction was 
first considered in the context of studying spectral statistics in [22] and generalized in [28] 
to directed graphs. On each bond the waves travel in both directions. Therefore, at any 
given moment, the system is fully specified by a 2S-dimensional vector of wave amplitudes 
a, indexed again by the directed bonds. Scattering at vertices is described by a unitary 
matrix T, having the property that 

^(fl,«2)(f3,f4) =0 if ""2 / ^3- 

The free propagation results in the amplitude (here 6 is a directed bond) being multiplied 
by the phase factor e'''*^''. Altogether we arrive to the quantum evolution operator U{X) = 

Thus, in both constructions one ends up with a unitary matrix ^(A) = e''^^5'o. This 
matrix specifies the eigenvalues {A„} via the equation 

det[I-e'^^So] = 0. (6) 



Actually, it was noted in [29] (see also [11[27]) that the multiplicity of A = as a root of 
([6]) can be different to its multiplicity as an eigenvalue of the Laplace operator, but that for 
all positive eigenvalues, the multiplicities coincide. To avoid this ambiguity we will omit 
the zero eigenvalues from the spectra that we consider, i.e. in our counting, Ai will be the 
first strictly positive eigenvalue. Since the multiplicity of the zero eigenvalue is necessarily 
finite (see corollary 14.21 below) this omission will not affect any spectral statistics. 



2.3 Spectral quantities for quantum graphs 

In various sources the notion of the "spectrum cr(G) of the graph G" can refer to:- 
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1. the eigenproblem ([T]) and thus solutions {A„} of ([6]), 

2. the eigenphases of the matrix S{X) for an arbitrary A, i.e. to the set of 2B numbers 
{9j} such that e'^^ is the j-th eigenvalue of S{X). 

To distinguish the two notions of the spectrum we will refer to the first definition above 
as the A-spectrum, and to the second as the 0-spectrum. 

In a similar way, the "eigenvector" of G can refer to one of three objects:- 

1. the function u{x) that solves ([1]), subject to boundary conditions, for some A^ in the 
A-spectrum, 

2. the eigenvector of 5'(A„) corresponding to the eigenvalue 1, i.e. to the solution c of 
S{Xn)c = c. These eigenvectors will be denoted by 0^, 

3. any of the 2B eigenvectors of S{X) for arbitrary A, denoted by tpj{X), j = 1, . . . , 2B. 

There correspondence between first two types of eigenvectors is given by formula ([3]). A 
heuristic formula which connects properties of the second and the third types of eigen- 
vectors was suggested in [lOl Eq. (5)]. In this article we give a rigorous derivation of the 
formula. 

3 Main results 

Our results concern equivalences between statistical properties of two types of eigenval- 
ues and eigenvectors of quantum graphs. These results are derived as an application of 
Proposition 14.41 introduced in Section [H 

3.1 Eigenvalue statistics 

Starting with the seminal work of Kottos and Smilansky [HI [22], quantum graphs became 
a popular model of quantum chaos. One of the more pertinent questions of quantum 
chaos is the universality of the eigenvalue correlations among systems of certain type. 
On graphs, [21j showed some preliminary numerical evidence that eigenvalue spectrum of 
graphs follow the general prediction [30^ i31j which says that the spectrum of classically 
chaotic system should have correlations typical associated to eigenvalues of large random 
matrices. Persistent deviations from the predicted behavior were found in star graphs \12\ 
l32] and Tanner [7] proposed a precise condition on the graphs to follow the random matrix 
theory prediction. This question was then attacked analytically by various methods, with 
results reported, in particular, in [U [5l [9]. For more information we refer the reader to a 
recent review [1]. 

While the original interest was in level statistics, i.e. the statistical functions of the 
spectrum {A„}, most of the analytical studies were in fact concentrating on the eigenphase 
statistics of the 0-spectrum. There are heuristic reasons why the corresponding statistics 
should coincide, but the formal link between the two has hitherto not been explored. 

In this article we prove that, in the correct limit, the two statistics are equivalent. 
We take the example of nearest-neighbor spacing distribution, i.e. the distribution of gaps 
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{Xn — A„_i} on one hand and the values of the functions {6j{\) — 0j-i(A)} on the other. 
Formahy, the level spacing distribution is given by 

1 ^ 

P(s)= hm ^5{s-iXn-Xn.i)), 

iV — >oc 1\ ^ — ' 

n=l 

where S is the Dirac delta-function. Mathematically, the distributions are defined via a 
family of test functions. 

Theorem 3.1. Let L be linearly independent over Q and L = (Li + . . . + Lb)/B denote 
the mean bond length. If h be a continuous function, then the limits 

1 ^ 

Px[h] := hm - Y^h (L(A„ - A„_i)) , (7) 

TV — >oo I\ ^ — ' 

n=l 

and 

Pe[h]=^mjf -L^/,(0^.(A)-^^._i(A))dA (8) 
-JO j^i 

exist and define bounded linear functionals of h. If we take the limit L — > ^o(l) !> • • • > 1) for 
some io > 0, while keeping L linearly independent over Q, the two functionals coincide. 
Namely, 

lim Px\h] = lim Pg\h], (9) 
where AL := Lmax — -^^min is the spread of the distribution of the bond lengths. 

Remark 3.2. If L is not linearly independent over Q, the limits ([7]) and ([8]) still exist 
but are highly sensitive to the changes in individual bond length. Thus the limits in 
equation ([9j) are, in general, singular. In particular when all bond lengths are equal 
(AL = 0), the A-spectrum is periodic and spacings in ^-spectrum are independent of A. 
The nearest-neighbor distribution in this case is highly degenerate. 

Remark 3.3. It is straightforward to extend the result to r-th nearest neighbor spacing 
distributions, i.e. the distributions of A„ — Xn-r and 6j{X) — 9j-r{X). Moreover, when all 
r-th nearest neighbor spacing distributions coincide, so do other statistical functions such 
as the n-point correlation functions. (See, for example, equation (6.1.39) of [33].) 



3.2 Eigenvector statistics 

The equivalence between the statistics of the A-spectrum and the 0-spectrum can be 
extended to eigenfunction statistics. 

To proceed, we need to introduce some notation in addition to that introduced previ- 
ously in section [2^31 By •0j(A) we will denote the j-th eigenvector of 5(A). As before, L 

denotes the vector of the bond lengths. The total length of the graph is .if = X^^^ Lf, 
and the average bond length, 5£ jB, is denoted by L. 

Let ^ be a X 25 matrix ("observable"). We denote An = {(fiJAlcf)^) and Aj{X) = 
{tpj{X)\A\ipj{X)), which corresponds to the "expected value of the observable A". We 
define the bond-length observable by L := diag(L,L) and set L{Xn) = {4>„\L\4>„). 
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Further, let D denote the 2B x 2B random unitary matrix 

Z):=e'^, 

where X is a random vector distributed uniformly on := [0,27r)^, and Aj{D) = 
{iJ}j{D)\A\ipj{D)) be the expectation of ^ with respect to the j-th eigenvalue of the matrix 
DSq. By E^(-) we denote averages with respect to the distribution of D. 

One of the central results in the study of statistical properties of eigenvectors is "quan- 
tum ergodicity" . Quantum ergodicity is the property of almost all eigenvectors to equidis- 
tribute [34 1 135 ^ [36], In our context, this corresponds to vanishing of the variance of either 
An or Aj(X) in some limit. We shall not discuss here the important question of which limit 
is appropriate; we refer the interested reader to [10]. Instead we prove that the variances 
(as indeed all other moments) of An and Aj{X) are equivalent. 

Theorem 3.4. Let L be linearly independent over Q. Then, for all m > 0, 

N g I(a37I = S X 21 g ^^<^'"' = ^° (21 g •*^<°''" j • '1°' 

The weighting on the left-hand side of (|10|) shows that there is an exact equivalence 
of moments only in the (additional) limit AL 0, as was taken in [32\ 137] . However, 
provided AL is bounded away from zero and infinity, the centered moments would ap- 
proach the zero limit either simultaneously or not at all. Another important feature of 
equation ([TO]) is that its second part relates the moments of the eigenvectors of a fixed 
graph to the averaged properties of an ensemble of random matrices. While such averaging 
was not necessary in the proof of quantum ergodicity of models considered in [10], it is 
expected to be helpful in more general circumstances. 

4 Main tools 

4.1 Graph spectra and torus flow 

Our methods are based on an idea going back to Barra and Gaspard [38], who suggested 
to view the eigenvalues An as the Poincare return times of a flow to a hypersurface defined 
by an extension of equation Q . 

Let T"^ be the torus of side 2tt in B dimensions. Define the surface S C T-^ by 

S := {x : detie'^'S'o - /] = O} , (11) 

where we use notation ([5]) for exponential of a vector. 
Further, define a flow (/>((xo) on T-^ by 

4>t{^o) = xq + tL mod 2tt. (12) 

If the components of L are linearly independent over Q then the flow is equidistributing 
on the torus. 

Since (pxW — mod 27r, equation ([6]) can be rewritten as 

det[e''^^(°)5o - /] = 0, 
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Figure 1: A cartoon of the flow piercing the surface E. The direction of increasing k is 
indicated by an arrow, and eigenvalues are indicated by Httle circles. 



or simply as (j)x{0) £ S. Thus the times t = An, n G Z of intersections of the flow (pt{f)) 
with the surface S give the set of points in the quantum spectrum of a graph (see figure 

The idea of Barra and Gaspard [38] was to use the ergodicity of the flow ([12]) to 
compute averages like ([7]) as an integral over the surface S. In this paper we put their 
idea on a mathematical footing (in particular, tackling the problem of using the ergodic 
theorem to integrate over S which is measure zero subset of T^) and extend it to apply 
to the problems described in Section [S] 



4.2 More about intersections with S 



The following result goes back to at least ^22^ Eq. (70)]. For completeness, we provide it 
with a proof. 

Lemma 4.1. Let U be a 2B x 2B unitary matrix, and let D{t) := e"^*^'"') for fixed xq. 
Let e'^^*^ be an eigenvalue of the unitary matrix D{t)U and u(t) be the corresponding 
normalised eigenvector. Then, 

c\f) 

- = (u(t)|L|u(i)), (13) 
where L = diag(L, L). In particular, 

Lmin ^ ^ ^max- (14) 

Proof. Since D[t)U is unitary and analytic, the function 9{t) can be chosen real analytic, 
and so there are no problems differentiating the eigenvalue equation, 

D(t)C/u(i) = e'^Wu(t), (15) 

with respect to t to obtain 

iLD{t)Uu{t) + D{t)Uu'{t) = i0'(Oe'^(*V(t) + e'^Wu'(t). (16) 
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Multiplying by u(t)"^ (the conjugate transpose of u(t)), 

ie'^(*)u(t)tLu(t) + vL{t)^D{t)Uu'{t) = i0'(t)e^^W, (17) 

where we used that u(t)^u(t) = 1 and u(t)"fu'(t) = 0. To handle the second term on the 
left-hand side we notice that 

u{t)^D{t)U = ({D{t)U)^u{t)y = fe-'^Wu(t))^ = e'^Wu(t)t (18) 



and so 

u{t)^D{t)Uu'{t) = e'^Wu(t)tu'(t) = 0. 



Thus, 

dt 



r\f) ^ 

u(t)tLu(t) = Y,HHt),\' + \n{t)j+B\% (19) 



leading to the desired bounds. □ 
We can fix bounds on the rate of crossing of the flow with E based on the fact that 

det[e^'^*(^o)5o-/] =0 

if and only if one of the eigenphases 6{t) is equal to a multiple of 2tt. Since all 9{t) 
are increasing functions and Lemma 14.11 gives us an estimate of the increase, we get the 
following corollary. 

Corollary 4.2. If t varies in an interval of length — then the flow </<t(xo) intersects 

the surface S at least 2B times. 

2tt 

If t varies in an interval of length — then the flow (j)t{^o) intersects the surface S 

at most 2B times. 

2-K 

On average, the flow (/)j(xo) intersects the surface 2B times in an interval of length 
_ L 

where L is the average bond length. This is a consequence of the asymptotic density of 
the A-spectrum, the Weyl's law for graphs. This classical result has been known on graphs 
for some time. A particularly elegant proof has appeared recently in j39l Theorem 1]. 

Proposition 4.3. Let N{A) := ^^{0 < An < A} where spectral points are counted with 
multiplicity. Then 

N{A) ~ —A as A ^oo, (20) 

TT 

where =Sf := Li + ■■■-)- L^- 

We remark that, in fact, the stronger asymptotic 

iV(A) = :^A + 0(1). 
vr 



is proved in [39], however we shall only make use of (|20j) . 
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Figure 2: The region in which the "thickened" function $e is supported. If two branches of 
S come close together, the function takes the sums of the values on the different branches. 
In the figure this region is shaded in a darker tone. 

4.3 The main proposition 

Let ^' be a bounded function, defined on S. We extend it to the whole of T^, by the 
following "thickening". Let e > and define 

ci.,(x) ■.= ^m^ {4=x+tL: |t|<e/2}- 

(21) 

Thus takes the value ^{$,) at all points within a distance at most e/2 from ^ € S in the 
direction of the fiow. Should two or more different branches of the surface S come close 
together, or intersect, the function takes the sum of the values of $(^) on the branches 
(see figure [2]). The sum in (|2ip is finite as a consequence of corollarv l4.2[ 

We permit the function $ to be undefined on a set U CTi, provided that U is small in 
the following sense. The projection of the set U onto the B — 1 dimensional hyperplane 
P := {x S : xi = 0} is required to have measure zero: 

meas {y £ P : 3t, < t < 27r/Li, 0t(y) £U} = 0. (22) 

This ensures that the discontinuities of $e have measure zero with respect to S. 
Our main tool is the following proposition. 

Proposition 4.4. Let L be linearly independent over Q and let <l> be defined on all of 
T, except possibly for a set U satisfying ()22p . and be positive and bounded. Denote by 
tn = tni^o) the consecutive intersection times of the flow </>t(xo) with the surface S. Then 
for any xq € and e > 

hm - V$((/>i„(xo)) = ^-— ^ / ci,^(x)dx, (23) 



7V-.ooAf^ ^^^"^ d(2n)B^J B 

n=l ^ ' 

where d := lim„-^oo n/tn = ^ /it is the mean density of the intersections t^- 

Remark 4.5. Since we want the result for every point xq € T'^, simple ergodicity of 
the fiow (pt would not be enough. Instead we shall use the Weyl's equidistribution which 
applies to intervals (and Riemann integrable functions) rather than Lebesgue measurable 
sets (and integrable function). Namely we use the following lemma, which follows directly 
from the results of Weyl [40] . 



9 



Lemma 4.6. Let G : — > M 6e a Riemann integrable function. Then, for every xq G T^, 

lim ^ r GiM^o))dt = -^ [ G(x)dx. (24) 

If the integral on the left-hand side does not exist, equation (^) should he understood in 
terms of upper and lower Darboux integrals (limits of Darboux sums), 

lim lL.D/^G(<^t(xo))dt= lim iu.D.rG(0t(xo))dt = -^ / G(x) dx. (25) 

Proof. The proof follows the standard procedure, see for example [4H Theorem 1.1]. Let 

= (J Rj (26) 



be a finite partition of T into well-behaved sets, such as rectangles. Then 



G(x)<G(x) :=^Iij,(x) sup G(y), 



where Ia is the indicator function of set A. By [40\ Satz 5], 

1 \R 



1 f 

lim - / lR.(6t(^o))dt 



the volume of the set Rj. Therefore, 



limsup;^U.D./ G((/.t(xo))dt < hm i / ^((/.^(xo)) dt = max G(y), 

which is the upper Darboux sum for the integral |T|~^ J^g G(x)dx. Since G is Riemann 
integrable, the Darboux sums converge to the value of the integral in the limit as the 
partition (|26p becomes finer. Analogous lower estimates lead to (j25p . □ 



Proof of Proposition \4-4\ Since the function is bounded and the discontinuity set has 
measure zero, <l>e is Riemann integrable and we can apply Lemma l4. 61 According to the 
way that $e is a thickening of width e of $ on S, we have that 

- $(</>, Jxo)) = —j^ ^,{M^o))dt + 0^-j , (27) 

as ^ oo. The error term comes from a possible overlap of the integration range with 
crossings other than tn, 1 ^ ^ A^- Corollary 14.21 guarantees that the number of such 
crossings is uniformly bounded. 

Now we let A^ ^ oo which entails T := t^ + e/2 ^ oo and apply Lemma 14.61 

hm l^<l>(<^,Jxo))= hm ^ hm 1 / '^M^o))dt (28) 

N-^co 1\ ^ — ' N^oo iVe T—*oo 1 In. 

n=l 

d(27r)^e Jjb 

Curiously, equation (j28|) provides a proof that the limit of n/t„ exists and is independent 
of Xq. Indeed, set <l*(x) = 1 and observe that two other limits in equation (j28p obviously 
exist. Thus, the limit of the sequence (n/t„) is the same as for the case xq = 0, which is 
covered by proposition 14.31 giving the limit ^/vr. □ 
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5 Applications 



With the tools developed in Section H] we can prove Theorems 13.11 and 13.41 While the 
proofs follow the same set of basic ideas, the eigenvector statistics proof is slightly simpler 
and thus we present it first. 



5.1 Eigenvector statistics 

In this section we will prove Theorem 13. 4[ We will do it by introducing a family of spectra, 
{Xa,n} indexed by a G [0,2-ir). We remind the reader that were defined as solutions 
of the equation det[e'^'^5o — /] = 0. We extend this definition and denote by Xa,n the 
solutions of the equation 

det[e-'"e'^^So - /] = 0. (29) 

Thus, Aq^„L — q; G S, where a = a(l, . . . , 1) and S is the surface defined by (llip . In other 
words, Ao,n are the intersection times of the flow 0f(— ct) with S. 

Similarly to 0„ we can now define (p^ „ to be the eigenvector of eigenvalue 1 of the 
matrix e'"e'-^"'"^5o. Obviously, An = Xo,n and 0„ = 0o,n- We let Aa,n ■= (0a,nl^l0a,n)- 

Lemma 5.1. Let G be a bounded continuous and possibly nonlinear functional on C^^. 
For X E S define <I*(x) := G {4>y^, where 4>^ satisfies 

e'^'S'o^x = <^x- 

Then 

1 ^ 1 r 

]vg^(^-) = (2^X.'^^(-)^- 

Proof. First we remark that '&(x) is well-defined and continuous on S, except possibly at 
degenerate points of S. 

By definition, 4>an — 4>:x. with x = Ao^„L — ot. Therefore, 

G {(l}a,n) = ^ (Aa,nL - q) = $ {^x^^ J-Oc)) . 

Now we can apply Proposition 14.41 to conclude the proof. □ 



The important consequence of Lemma [5. II is that the moments of Aa^n are independent 
of a. 

Proof of Theorem [gT^l Let {e'^-''^'^)}^f be the eigenvalues of the matrix e'^^^o, with 9j{X) 
chosen to be real analytic. As a consequence of lemma 14.11 9j has a smooth inverse. 

Using a 27r-periodized Dirac delta function, 62tt we can trivially write 

AjiXr = r hn{a - e,{X)) A,{Xr da. 
Jo 
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Then 

rA 2B r2n 2B M 



Jo Jo j=iJo 

=/. SI™, »'w-'(f)) 



0<A„.„<A 



We have chosen jQ^„ so that Oj^^„{Xa,n) = o mod 27r. Therefore t/jj^ ^ (Aa,n) is what we 
denoted by 4>a,n- Further, by Lemma Wl\ 9j{Xa,n) = {(t>a,n\^4>a,n) ■ We set 



and obtain 



^0 j=i Jo A..„<A 



(32) 



Now we divide by 2BA and take the Hmit A ^ oo. By proposition 14.31 A^(A)/A J^'/tt 
in this hmit. By the dominated convergence theorem we can take the hmit inside the 
a-integral to get 

lini - / — y 4.(A)'»dA = ^ / hm —-- V G(0„J da 
A^^AJo 2^4- ^ 2B7tJo A^ooiV(A)^^^ ^^"'"^ 



JV(A) 



1 A A!" 



= hm — y (33) 
^^°°^;^^(An)/i^ 

where we used the fact that, by Lemma [5. H the hmit inside the integral is independent of 
a. 

To prove the second equahty in (fTOl) we apply Weyl's equidistribution theorem to the 
limit on the left in (j33p . This is justified since eigenvector components of e^^So vary 
continuously in x except possibly at the points where the eigenvalues of e'^So are non- 
simple. The set of such points is a C-analytic set and will have either measure or full 
measure (see the discussion in Section 15. 3p . The latter situation is incompatible with 
condition (1221). □ 
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5.2 Spectral statistics 

We will use the previous results to prove that the spectral statistics of the matrix e'^'S'o 
averaged over converge to the empirical spectral statistics as the interval from which 
bond lengths are drawn shrinks in size. 

5.2.1 A few more definitions 

Let d : R+ be the "time" of the next crossing of the surface S, 

(i(x) := inf{t > : (/>t(x) G S}. (34) 

For X e T-^ let the eigenvalues of e'^5o be e'^^'^^^ . . . , e^^^^^^\ ordered so that 

< ^2b(x) < ^2B-i(x) < • • • < ^"2(x) < ^i(x) < 2tt. (35) 

This ordering ensures that if x € S then ^i(x) = 2it. We use the hat in notation 6j to 
distinguish it from the eigenphases 9j{X) of the matrix e'^^5o which can be chosen real 

analytic with respect to A. Of course, the sets {%(A)} and {(f)x{0))^ coincide modulo 
27r. 

Define the eigenphase spacing functions by 

(71 (x) :=^"i(x)-^2(x) 

0-2B-l(x) := 02B~l{x) - 023 (x) 

C72ij(x) := e2B{x) +2tt- ^i(x). 

These definitions are illustrated in figure [3l 

We will use shorthand crj(A) for aj{(px{0)) and d{X) for d{(p\{0)). 

5.2.2 Equivalence of spacing distributions 

Proof of Theorem \3.1[ Let hhe a continuous test function. We want to compare 

A 2B N 

lim - / — V/i(fT,(A)) dA and lim — V /i (Ld(A„)) . 

■^^ j=l n=l 

Repeating the arguments of equations (|31|) and (|32|) we obtain 

rV/i((j,(A)) dA= VM«-^r-(A))/i(cT,„(A)) dadA 

Jo Jo Jo 

= / / d{X-Xa,n)—^, — TY^dAda, 
•''0 Jo ^j.J^> 

= E ^''('Aa..„(-«)) da, 
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o-i(0a(x)) 



O-2(0a(x)) 



O-3(0a(x)) 



A = Ai A 

Figure 3: The eigenangles of e^'^^^^^So and notation of Section [5.2.11 



where, for x G S, we have defined 

$'^(x) 



dt 



(36) 



t=o 



where the derivative is taken from the left, since 9i is discontinuous on S. The function 
$'^(x) is weh defined if at x we have 62 < Oi = 27r, i.e. the eigenvalue 1 of e'^S'o is simple. 
Now we divide by 2BK, take the limit A — > 00 and use Proposition 14.41 to get, analogously 
to (IMD, 



On the other hand, again by Proposition 14.41 we have 



«>^(x) dx. 



1 ^ 1 /■ 

Jnn^-X;/.(Zd(A„))=^y^^cl>f(x)d> 



where 



$'='(x) = /i (Ld(x)) . 



(37) 



Now we use that L^>'^(x) — <I>'^(x) ^ as AL 0. This is proved in lemma [5^2] below, 
implying 

I 



[ $-(x)dx--l^ / cl>f(x)dx^0. 



(38) 



Since $'^(x) is clearly continuous in L, the individual limits exist too, concluding the 
proof. □ 
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Lemma 5.2. Let h be continuous. Then, for x E S we have 

^ < d(x) < (39) 

and, therefore, 

Z^>'"(x) - $'^(x) ^ (40) 
uniformly on T, as AL = L^ax — -^^min 0. 

Proof. Estimate (|39p follows from lemma 14.11 Indeed, (i(x) is the time required by 
01 ((jit (x)) to reach 27r from its initial value of 27r — cji + at the time t = +0, see Figure [3l 
and lemma HTTl provides an estimate on the derivative of 9i. 

We can now write ()39p as 

ai(x) (j^ -l) < Ld{^) - ai(x) < ai(x) (-^ - l) . 

Since both d and ai are bounded, we have Ld{x) — (Ti(x) uniformly in x as the 
spread of the lengths decreases. The function h is continuous on a compact set, therefore 
it is uniformly continuous and h{ai{x)) — h{Ld{x)) — > uniformly in x. Finally, the 
denominator in the definition of ^I^*^, see equation is bounded by -Lmin and -^max by 
lemma WA\ and thus converges to L. All of the above put together imply ()40p . □ 



5.3 Discontinuities of and ^'^ 

When defining the function ^ in Section 15.11 and the functions and ^I''^ in Section 15.21 
we implicitly assumed that they satisfy condition ()22p on the set of their discontinuities. 
In this section we explain why can we expect it in general and how to rectify the situation 
when the definitions of function $ produce too many discontinuities. 

The surface S is a C-analytic set, see [42 ^ Definition 6, Chapter V]. Thus we can 
apply |32i Proposition 18, Chapter V], to it, which splits S into an analytic manifold 
of dimension B — 1 and the remainder U. The remainder has a smaller dimension and 
obviously satisfies condition ([22]) . 

Consider a connected piece of the {B — l)-dimensional manifold. If on this piece the 
zero of the defining equation 

det[e'<^*(^o)5o-/] =0 

is simple, there is a unique eigenvector 0^ (see Section [5T]) . continuous with respect to x 
on the manifold. The function $ is well defined. 

The situation when on all {B — l)-dimensional manifolds the eigenvalues are simple 
is expected to be generic. However, a proof exists only for a special case of graphs with 
Kirchhoff matching conditions [43]. Proving the simplicity (in a generic sense) of the 
solutions to 

det[e'''5o - /] = 

for general 5*0 thus remains an important open question. To complicate the picture, there 
are counter-examples with persistent eigenvalues of multiplicity two: these are provided 
by graphs with a looping bond separated from the rest of the graph. Moreover, there are 
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indications that the first-order perturbation used in [43j will not produce the result for 
general Sq. 

Fortunately, when the multiplicity is greater than one on a (i? — l)-dimensional sub- 
manifold of S, we can still make the functions ^ continuous. First we note that the 
multiplicity must be constant on the whole manifold. For simplicity we assume the mul- 
tiplicity is two. The eigenspace of the eigenvalue 1 of the matrix e^'^^^^^^So is continuous 
and we can choose a continuous basis i ™d 0x2- The function $ of Section [5TT] can 
be then defined as 

on the problematic manifold. 

The function of Section [52] can be simply defined as h{Ld{ii.)) + h{0). To define 
function in a continuous and meaningful manner we note that if ^i(x) = 02(x) = 27r, 
then the phases also coincided before the flow hit S, i.e. for small negative t, 0i(0f(x)) = 
02('/'t(x)). Thus the derivative in the denominator of equation (j36p does not depend on 
which 6 we take and one can set 

^^^^^ ^ /i(ai(x)) + MO) 
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